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Abstract 


本 文 介 绍 了 环形 加 速 器 中 的 粒子 动力 学 ， 环 形 加 速 器 中 带电 粒子 的 运动 分 为 横 
向 运动 与 纵向 运动 。 本 文 的 题材 来 自 于 加 速 器 的 理论 与 设计 ， 我 们 试图 将 物理 问题 
与 数学 方法 有 机 地 结合 起 来 讨论 环形 加 速 嚣 中 的 粒子 动力 学 。 一 方面 ， 物 理 问题 的 
提出 推动 了 数学 的 发 展 ， 男 一 方面 ， 数 学 的 发 展 使 得 物理 问题 可 采用 更 加 高 级 的 数 
学 工具 来 简化 物理 问题 的 解决 ， 从 而 解决 物理 问题 更 加 容易 。 本 文 的 特色 是 运用 哈 
密 顿 力学 处 理 环形 加 速 器 中 的 粒子 的 横向 运动 和 纵向 运动 。 
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1 ”引言 


环形 加 速 器 的 概念 是 上 个 世纪 五 十 年 代 提 出 的 ，1952 年 E.Courant,H.Synder,and 
S.Livingston 提 出 了 交 变 梯度 聚焦 的 概念 ， 交 变 梯度 聚焦 也 就 是 强 聚 焦 原理 。 强 聚焦 原 
理解 决 了 回旋 加 速 器 需要 大 型 磁铁 的 困难 ,回旋 加 速 器 向 更 高 能 量 发 展 需要 的 磁铁 越 来 
越 大 ， 它 是 弱 聚 焦 加 速 器 。 加 速 器 中 束 流 靠 四 极 磁 铁 完 成 ， 单 块 的 四 极 铁 在 一 个 横 癌 
方向 聚焦 ， 在 另 一 个 横向 方向 散 焦 ， 单 块 的 四 极 铁 无 法 同时 在 两 个 横向 方向 聚焦 ， 但 
是 四 极 磁 铁 的 组 合 在 整体 上 可 以 完成 两 个 横向 方向 的 聚焦 。 假 定 两 块 符号 相反 的 四 极 
磁铁 的 等 效 焦距 为 方 , 户 ,这 两 块 四 极 磁铁 间隔 距离 4 放置 ， 那 么 其 等 价 的 四 极 磁铁 的 焦 
BEF 
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fi fo Ss Ad RATT DARE > 0。 那 么 这 两 块 间隔 放置 四 极 磁铁 整体 上 在 两 个 横 
向 方向 都 是 聚焦 ， 这 就 是 强 聚 焦 原 理 ， 亦 即 交 变 梯度 聚焦 原理 。 环 形 加 速 器 在 五 十 
年 代 、 六 十 年 代 获 得 了 发 展 ， 它 是 加 速 器 发 展 史 上 的 一 个 飞跃 。 在 直线 加 速 器 中 ， 带 
电 粒子 的 纵向 运动 稳定 性 与 横向 运动 稳定 性 不 兼容 ， 纵 向 聚 束 和 横向 聚焦 无 法 同时 完 
成 ， 所 以 也 一 样 采取 强 聚 焦 原理 对 带电 粒子 束 进行 聚焦 ， 也 就 是 用 四 极 铁 的 组 合 进行 
横向 聚焦 。 如 果 没 有 强 聚 焦 原理 ， 直 线 加 速 器 不 可 能 把 带电 粒子 束 加 速 到 高 能 量 。 而 
在 环形 加 速 器 中 ， 带 电 粒子 的 加 速 是 在 高 频 腔 中 补充 或 加 速 能 量 ， 带 电 粒 子 束 在 环形 
周期 磁 聚 焦 结构 中 依靠 强 聚焦 原理 得 到 了 很 好 的 横向 聚焦 ， 带 电 粒 子 束 才能 反复 的 通 
过 高 频 腔 得 到 补充 或 加 速 能 量 。 


环形 加 速 器 主要 由 二 极 铁 、 四 极 铁 、 六 极 铁 、 高 频 腔 所 组 成 ， 带 电 粒 子 在 二 极 
铁 中 偏转 ， 在 四 极 铁 、 六 极 铁 中 聚焦 ， 在 高 频 腔 中 加 速 ， 带 电 粒 子 在 坏 中 运动 、 储 
存 。 环 形 加 速 器 主要 有 储存 环 (storage ring)、 增 强 器 (booster)、 阻 尼 环 ( damping ring 
)、RLA( Recirculating linear acelerator ) 环 /ERL( Energy recovery linacs ) 环 -返回 型 直 
线 加 速 器 返回 坏 。 在 储存 坏 中 带电 粒子 作 环 形 运 动 ， 储 存在 环 中 ， 同 时 在 高 频 腔 中 得 
到 加 速 补充 由 于 辐射 而 损失 的 能 量 。 增 强 器 在 储存 带电 粒子 的 同时 在 高 频 腔 中 反复 加 
速 带电 粒子 ， 使 带电 粒子 从 低能 加 速 到 高 能 。 阻 尼 环 在 储存 带电 粒子 的 同时 通过 扭 摆 
磁铁 阻尼 带电 粒子 束 的 发 射 度 。RLA 环 /ERL 环 在 储存 带电 粒子 返回 的 同时 通过 波 荡 器 
兼作 光源 。 

储存 坏 是 最 主要 的 环形 加 速 器 ， 主 要 应 用 于 高 能 对 撞 机 与 同步 辐射 光源 。 用 正 
负电 子 作 对 撞 实 验 研 究 粒 子 物理 是 20 世 纪 50 年 代 提 出 的 ，20 世 纪 60 年 代 建 成 了 最 初 
的 儿 台 正 负 电子 对 撞 机 ， 成 功 进 行 了 正 负 电子 对 撞 。 它 们 有 Frascati 的 250MeV 储 存 
环 AdA(1962 年 ), 斯 坦 福 的 500MeV 储 存 环 CBX (1963 年 ) 前 苏联 新 西伯 利 亚 的 140Mev 储 
存 环 VEP-1(1963 年 )。 同 步 辐射 光源 是 在 正 负 电子 对 撞 机 的 基础 上 发 展 起 来 的 ， 大 多 数 
对 撞 机 兼作 同步 辐射 ，20 世 纪 70 年 代 以 后 开始 建造 专门 的 同步 辐射 光源 。 同 步 辐射 装 
置 经 历 了 一 代 、 二 代 、 三 代 、 四 代 , 是 非 相 干 辐射 光源 . 一 代 同 步 辐射 装置 是 寄生 在 高 
能 物理 装置 上 兼用 的 同步 辐射 装置 ( 20 世纪 70 FRAR). 二 代 同 步 辐 射 装置 是 大 量 放 
入 搬入 件 的 专用 光源 ( 20 世纪 80 年 代 以 来 ). 而 三 代 同 步 辐 射 装置 是 优化 的 , 低 发 射 度 的 
专用 的 同步 辐射 光源 ( 20 世纪 90 年 代 以 来 ). 而 目前 世界 上 普遍 研究 、 设 计 、 建 造 、 运 
维 的 是 第 四 代 同 步 辐射 光源 -衍射 极限 环 ， 它 的 发 射 度 、 亮 度 都 达到 了 极限 。 


2 ”哈密 顿 力学 


1788 年 拉 格 朗 日 写 了 著作 《分 析 力 学 》， 完 全 用 数学 分 析 的 方法 来 解决 力学 问题 ， 
无 须 借助 以 往常 用 的 几何 方法 ， 并 逐步 发 展 成 为 一 系列 处 理 力 学 问题 的 新 方法 ， 即 所 
请 的 分 析 力 学 。 拉 格 朗 日 是 用 n 个 独立 变量 来 描写 力学 体系 的 运动 ， 所 以 和 牛顿 运动 方 
程 一 样 ， 是 二 阶 常 微分 方程 ， 我 们 把 这 组 方程 称 之 为 拉 格 朗 日 方程 。1834 年 哈密 顿 提 
出 ， 如 果 用 坐标 和 动量 为 独立 变量 ， 那 么 运动 方程 (微分 方程 组 ) 有 二 阶 降 为 一 阶 ， 
但 方程 式 的 数目 翻 了 一 倍 ， 变 为 2n 个 方程 式 ， 这 组 方程 我 们 称 之 为 哈密 顿 正则 方程 。 
哈密 顿 正则 方程 在 理论 上 有 重要 意义 ， 它 是 从 经 典 力学 走向 量子 力学 的 桥梁 ， 在 理论 
物理 的 各 分 支 中 有 广泛 的 应 用 。 哈 密 顿 在 1843 年 又 提出 运用 变 分 法 的 哈密 顿 原理 来 描 
写 力 学 体系 的 运动 ， 它 和 和 牛顿 方程 、 拉 格 朗 日 方程 、 哈 密 顿 正则 方程 等 价 。 


2.1 拉 格 朗 日 力学 与 喻 密 顿 力学 

假定 具有 nn 个 自由 度 的 力学 体系 的 坐标 变量 gi,i = 1,.… ,n, 广 义 速 度 变 量 G; = fi = 
1,.… ,n, {为 时 间 变 量 ， 那 么 拉 格 朗 日 函数 L(g;, Gi,t) 是 2n 个 动力 学 变量 的 函数 。 运 用 力 
学 的 变 分 原理 (哈密 顿 原理 ), 则 有 变 分 方程 
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由 以 上 变 分 方程 可 得 拉 格 朗 日 方程 为 

d ‚ðL OL 

dag) Bq ©) 
引入 广义 动量 变量 p, = 22, EN BL WARMER (qi do ERBER (a po t), HELE 
密 顿 函数 为 


H(q,p,t = (gq,p,t), (4) 
哈密 顿 通 数 的 微分 为 
OH OH OH 
H= | i | 2) 
d dt dq OF, dp (5) 


同时 由 哈密 顿 函数 的 定义 ， 其 变 分 展开 为 


. OL OL .. OL 
dH = > Pidi + 2 Gidpi 一 i T = gA = pr (6) 
由 以 上 微分 形式 可 得 以 下 方程 组 
OH 
qi Op, (7) 
OH 


以 上 方程 组 即 为 哈密 顿 正 则 方程 组 ， 它 是 2n 个 一 阶 微分 方程 组 ， 同 时 有 红 = 中。 
带电 粒子 在 电磁 场 中 运动 时 ， 其 拉 格 朗 日 函数 为 


L= —moe"(1— 5)? — eb + eA 0, (9) 
Th mo, e 为 带电 粒子 的 质量 与 电荷 ，9, 4 为 电磁 场 的 电位 势 和 磁 矢 势 。 而 哈密 顿 函数 为 
H = eġ + c|(p — eA}? + me]?. (10) 


我 们 在 求解 哈密 顿 方程 时 ， 经 常 通过 一 些 正则 变换 ， 寻 找 运动 不 变量 ， 进 而 简化 哈 
密 顿 方程 。 正 则 变换 是 通过 系统 坐标 变换 来 进行 的 ， 对 于 坐标 变换 (gi, pi) 一 (Qi, Pi), 
其 构成 正则 变换 的 充分 必要 条 件 是 


bs Hil( (Qi, Pa t)] 这 RE H (qi, pi, t)] = dG, (11) 


qdG 是 一 个 完整 的 微分 。 对 于 正则 变换 ， 新 旧 坐 标 之 间 的 关系 为 


Qi = Qi(q 205 On Pls , Pn, t), (12) 
P; = P;(q, DO , Pnt). (13) 


变换 关系 由 生成 函数 G 求 得 ， 有 四 种 正则 变换 : 
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° Gi(q, Q, t) : Pi = i = 


s Gə(q, P,t) : pi = 3 Or 
° G3(Q, p, t): P,=-#,4=-%, 


e Gulp, Pt): qi = 一 se. Qi = $F, 
而 变换 后 的 哈密 顿 函数 万 ;为 


OG 


Hi(®,P,t) = H(qp,t) + a. (14) 


2.2 ”哈密 顿 力学 在 加 速 器 中 的 应 用 
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Figure 1. 环形 加 速 器 的 参考 坐标 系 


在 环形 加 速 器 中 ， 假 定 带 电 粒 子 的 参考 闭合 轨道 在 一 个 平面 上 ， 或 者 说 在 水 平面 ， 如 
上 图 。 引 入 Frenet-Serret 坐 标 系 ， 参 考 轨 WE, s 为 参考 轨道 上 参考 点 到 起 始点 
之 间 的 轨道 长 度 ， t 为 轨道 切 向 单位 向 量 ， 力 为 轨道 法 向 单位 向 量 ，? 为 轨道 平面 法 向 单 
位 向 量 ， 角 ,5,1 构成 一 个 右手 正 交 坐 标 系 。 那 么 围绕 参考 轨道 的 粒子 向 量 可 表达 为 


r(s) =ro(s) + zñ + yb (15) 


有 如 下 特性 
必 (16) 
m Z -zi (17) 
ara (18) 
De (19) 


上 式 中 ，p 为 参考 轨道 曲率 半径 。 
在 坐标 系 (x,y,s) 中 , 假定 带电 粒子 的 正则 动量 为 户 , 作 正 则 变换 ， 定 义 正则 变换 生成 


CP Gas SP one AA, (20) 
对 于 坐标 系 (z, Vs) 的 正则 动量 为 
aG i 
paea EE (21) 
ac f 
p=- = Ph, (22) 
aG : 
ag oP (23) 


AAS = 0, 也 就 是 G 不 含 时 间 变 量 ， 在 环形 加 速 器 中 ， 横 向 运动 只 考虑 在 磁场 中 的 运 
动 ， 假 定 磁 矢 势 为 4, 在 坐标 系 (x,y, s) 中 类 似 有 


Ay = (1- —)A-4, (24) 
p 

A, = A-h, (25) 

A, = A-6. (26) 


从 而 ， 带 电 粒 子 在 磁场 中 的 哈密 顿 函 数 为 


H(q,p,t) = cl(ps — eA)” + (py — eAy)” 4 +- m2c2]3. (27) 


在 环形 加 速 器 粒子 动力 学 中 ， 我 们 往往 以 s 为 时 间 变 量 代 替 时 间 变 量 +， 由 哈密 顿 变 分 原 
理 有 以 下 变 分 方程 


dr dy ， ds B 
f wo t Puh t Psa — H)dt =0, (28) 
上 式 可 变换 为 
X | Y = | = 
5 fips fe Py as H oe | ps)ds 0. (29) 


由 上 式 可 知 ， 新 的 正则 变量 对 为 (z,p,), (y,Py),(t,-H), 而 哈密 顿 函数 为 = —p,, 时 间 
变量 为 哈密 顿 函数 已 = 一 六 为 


(于 一 mnac ) 一 (pz—eAz) 一 (py —eAy)7] 2 
(30) 


Tl 
FG a sb Px, Py, —H, s) SPs = Saga a 


从 而 哈密 顿 方 程 组 为 


=Z, (81) 
y= (2) 
i= 下- (33) 
n=- (34) 
r=- (35) 
Ht = (36) 


式 中 人 代表 至。 在 环形 加 速 嚣 中， 带电 粒子 在 静 磁 场 中 运动 ， 带 电 粒 子 在 二 极 铁 、 四 极 
铁 、 六 极 铁 中 运动 ， 粒 子 运动 的 哈密 顿 函 数 瑟 不 是 时 间 的 显 函 数 ， 不 舍 时 间 变 量 t, 所 
以 五 是 运动 不 变量 ， 是 常数 。 所 以 哈密 顿 方程 组 由 三 维 位 形 空间 约 化 为 二 维 维 形 空间 。 


2.3 ”多 级 磁铁 的 展开 与 模 向 运动 方程 


环形 加 速 器 主要 有 各 种 磁铁 所 组 成 ， 分 立功 能 的 磁铁 得 到 了 广泛 的 应 用 。 环 形 加 
速 器 中 带电 粒子 的 横向 运动 主要 是 带电 粒子 在 各 种 磁铁 中 的 横 癌 运动 。 磁 铁 的 种 类 很 
多 ， 主 要 有 二 极 铁 、 四 极 铁 、 六 极 铁 、 八 极 铁 、 螺 旋 线 圈 等 。 二 极 铁 用 于 带电 粒子 的 
偏转 ， 使 它们 的 轨道 构成 一 个 闭合 环 。 四 极 铁 用 于 带电 粒子 的 聚焦 。 六 极 铁 用 于 色 品 
的 补偿 和 非 线性 聚焦 。 八 极 铁 用 于 Landau 阻 尼 ， 螺 旋 线 圈 用 于 抑制 电子 云 不 稳定 性 。 

加 速 器 磁铁 的 理论 与 设计 属于 静 磁 学 问题 。 由 于 加 速 器 的 轴 对 称 结构 特性 ， 磁 铁 的 
理论 设计 可 以 简化 为 二 维 静 磁 学 问题 。 忽 略 纵向 运动 方向 坐标 ， 也 就 是 粒子 运动 前 进 
方向 ， 假 定 三 维 坐 标 系 (zx,y), MATE AB AT FE, 我 们 有 


VBN (37) 
引入 磁 矢 量 势 4， 
B=VxA, (38) 
A RAAT IAI) A,. 由 静 磁 学 理论 , 我 们 有 
VxB=0. (39) 
由 上 式 可 引入 标量 势 V， 
se (40) 
在 坐标 系 (z,y) 中 , 磁场 分 量 B,, By 可 推导 为 
B=- = (41) 
=m- (42) 


假定 复 变 量 z = z + iy, 定义 复 变 函 数 4 


A(z) = As(a,y), +iV (2, y). (43) 
在 复 平 面 > A, V 满足 Cauchy-Riemann 条 件 , 所 以 复 变 函数 4A(z) 是 解析 函数 ， 可 做 以 
下 的 级 数 展开 
= ` Kaz (44) 
n=0 


由 复 变 函数 分 析 , 对 所 有 在 区 域 | z |< re 内 的 点 xz 收敛 。 在 园区 域 之 外 , 复 变 函数 不 再 是 
解析 的 。 收 敛 半径 可 看 作为 轴 心 到 磁铁 斩 或 线圈 最 近 的 距离 。 


Im(z) 


radius 


Figure 2. 复 平 面 收敛 区 域 


在 实际 应 用 中 , 我 们 可 以 先 求解 4A,,V, 进而 在 坐标 系 中 得 到 磁场 分 量 B,, By EEF 
RAAB AR, ABE Kn = An + iln WA 


As(z,Yy) + iV (z, y) = ome" = =$ AFi, (45) 
n=0 
从 而 可 得 
As = bız + ayy + bo(x? — y’) + azzy + ..., (46) 


AP, a AR, AACE RAE ERR, bz, az 代表 正四 极 铁 和 和 斜 四 极 铁 。 
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环形 加 速 器 也 可 以 近似 看 作为 轴 对 称 位 形 。 
0, As = 4,(zx,y)。 对 于 四 极 磁铁 ，A, 的 表达 式 可 写 为 


在 多 极 磁铁 的 展开 中 ，4, = A, = 


Do 及 
As =~, Se =) (47) 
eg ae eres ae age 
为 四 极 磁 铁 的 中 心 轴 ) , 哈密 顿 函 数 简 化 为 
K í 
(ep (48) 
哈密 顿 方程 组 为 
fea Ee 49 
a (49) 
E 50 
P (50) 
P, = poKz, (51) 
Py = PoKy. (52) 
由 哈密 顿 方程 组 可 得 横向 运动 方程 
x" — Kz=0, (53) 
y + Ky=0, (54) 
横向 运动 方程 就 是 常 说 的 Hi 方程 。 
3 Waz) 
环形 加 速 器 主要 由 周期 性 磁 聚 焦 磁 场 组 成 ， 由 上 节 ， 带 电 粒 子 的 横向 运动 方程 为 
2" + K,(s)0 = 0, Ka = 5 F Kals), (55) 
y” + K,(s)y = 0, Ky = +K (8). (56) 


式 中 /为 曲率 半径 ， 有 为 实际 聚焦 强度 。Kv(s), K 
下 ， 横 向 运动 可 以 统一 由 一 般 性 的 Hil 方 程 来 描述 


p Po” 


天 (s) 为 周期 函数 ， Sr 为 能 散 或 动量 离散 。 假 定 周期 为 L, K(s+L) = 
解 由 齐 次 方程 的 解 和 特 解 组 成 ， 其 解 为 


x(s) 


xv" + K(s = 


A 
= C(s)xo + S(s)2h + Dl, 
0 
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,(s) 为 周期 函数 。 在 考虑 能 散 的 情况 


(57) 


K(s). Hill 方 程 的 


(58) 


式 中 zo = zx(s0), 起 始点 为 80，C(s)zxo + 5(s)zb 为 齐 次 方程 的 解 ，D(s) 全 方程 的 特 
fo C(s), 5(s) 为 齐 次 方程 的 两 个 线性 无 关 的 解 ， 其 满足 


C" + K(s)C =0,8" + K(s)S =0. (59) 


Hill 方 程 解 中 的 D(s) 定 义 为 色散 函数 ， 其 满足 以 下 方程 


1 

D" + K(s)D = 一 一 . 60 
S K 

定义 Wronsian 函数 如 下 

G aS 
W =| C’ Ss! p (61) 
W 满 足以 下 方程 
“a 二 全 (62) 
5 s 


那么 W 沿 粒子 轨道 是 运动 常数 ， 选 择 在 起 始点 s0, Co = 1,Ch = 0; So = 0,5, = 1, M 
WW = 1。 色 散 函 数 的 解 为 


Dis) = sts) f Cds- C) f sds (63) 


x(s),2’(s), Avs) 与 起 始点 z(so)， x'(S0), Avo) > iil 的 关系 由 和 矩阵 代数 表达 如 下 


(2 ).-Ce so) CG) +2 Ub): (6 

(65) 
3.1 ”磁铁 的 传输 矩阵 
Cs D 


(E $).(6 S D ) menses. msee, Kosaman 
0 0 1 


断 常 数 ， 定 义 角 变 量 o = s/|K|, EK > 0,K < 0,K = 0 三 种 情形 下 ，Hil 方 程 的 
解 C(s), 5(s) 为 
COs —_sin 
L K > 0 eA, (G oe i veo 
CS —,/|K|sing cosp 
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上 式 可 写 为 


—/|K|sinhy cosh 


: C S 1 s 
— DE PSY = 
3. K = 0, HETI, ( œ s ) = ( 0 1 | 


h h 
2. K <0, 散 焦 情形 ， ( C s ) 2 | cosho sin A 


kon 
1. K >0, D(s = Sf} Sds—C fy Sds= zx (1 — cosy), Di(s) = a 
2. K <0, D(s) = = — coshy), D'(s) = yp 
3. K =0, D(s) = 0, D'(s) = 0. 
由 以 上 公式 ， 可 得 漂移 段 、 四 极 铁 、 二 极 铁 的 传输 矩阵 如 下 : 
漂移 节 1=0,k=0 > K = 0, 传输 矩阵 为 
1 ! 0 
M,=M,=[ 0 1 0 (66) 
0 0 1 
四 极 铁 > =0>K = k, WRK > 0, 那么 传输 矩阵 为 
cosh ——sinhy 0 
VIKI 
M:= | VERlsinhp coshp 0 (67) 
0 0 1 
cosy —s—sinp 0 
VIKI 
My = | —V|Klsinp cosp 0 (68) 
0 0 1 


四 极 铁 在 水 平方 向 散 焦 ， 在 垂直 方向 聚焦 ， 如 果 开 < 0, Mı, Mi 互相 交换 ， 四 极 铁 
在 水 平方 向 聚焦 ， 在 垂直 方向 散 焦 。 在 很 多 情况 下 ， 薄 透镜 近似 计算 能 简化 矩阵 
分 析 , 对 Lattice 分 析 很 有 用 。 四 极 铁 的 焦距 /了 = 十 远大 于 四 极 铁 的 长 度 f = 4 > 1, 
在 薄 透 镜 近 似 下 四 极 铁 的 传输 矩阵 为 

) (69) 


0 
0 
1 
0 
0 ) (70) 
1 


二 极 铁 k=0 > K = 5, RHUF, p= VIKIll = i 其 传输 矩阵 为 


p 


cosp psinp p(l -— cosy) 
M; = | 一 25 和 0 Cosy siny ) (71) 
0 0 1 
1 l0 
M, = | 0 1.0 (72) 
001 


DA _E Fe BAAS HR EY 1 a EE DRE AE EY TR AR FS E ST EE E SF FT EE FY HR 
只 ， 我 们 以 一 个 telescope 望 远 镜 系 统 (doublet) 作 算 阵 分 析 。 


水 平面 垂直 面 
Figure 3. 望 远 镜 系统 (doublet) 


望远镜 系统 是 一 个 FODO 磁 聚焦 结构 ， 广 泛 用 于 高 能 加 速 器 中 粒子 束 的 聚焦 和 成 像 , WL 
配 lattice 函 数 。 如 上 图 ， 轨 迹 平行 于 轴 心 的 粒子 束 ， 无 论 在 水 平方 向 还 是 垂直 方向 都 能 
得 到 聚焦 ， 因 为 它 的 聚焦 幅度 大 于 散 焦 幅度 。 在 注 透 镜 近 似 下 ， 假 定 聚 焦 透镜 和 散 焦 
透镜 的 焦距 为 ,一 了 , 望远镜 系统 的 水 平 传输 矩阵 为 


1 0 0 010 1 00 bas l 0 
M doublets = F 1 0 0 1 0 = 1 0 = E T 0 5 (73) 
0 0 1 0 0 1 0 01 0 0 1 


其 在 垂直 方向 的 传输 矩阵 为 
e a ee 
M aoublet,y = | = 1 ; 0 : (74) 
0 1 


无 论 在 水 平方 向 还 是 垂直 方向 ， 传 输 和 矩阵 元 素 WM2i 都 为 Mol = = 它 是 telescope 系 统 的 
折射 率 ，telescope 系 统 在 两 个 方向 都 是 聚焦 的 。telescope 系 统 的 实际 焦距 为 


2 


fdoublet = a (75) 
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当 粒 子 从 距离 telescope 系 统 a 的 地 方 入 射 ， 它 在 离 telescope 系 统 ! 的 地 方 成 像 汇聚 ， 那 么 
其 传输 矩阵 为 


“Mr =S = 0, M 可 写 为 


m 0 0 
| œ 2 0 , (78) 
0 0 1 
从 a 点 到 b 点 可 以 成 像 ，z = mag, m 为 放大 倍数 。 如 果 C’ = =a 一 0, 那么 它 是 一 个 望 
远 镜 系 统 telescope。Mi, My"? Mi = S = 0 的 条 件 是 
bP eRe 
Ma: q= eta (79) 
a +i- 
Mi pS ae ey (80) 


上 式 中 ，a,b 是 可 交换 的 ， 当 a 关 5 时 是 散光 情形 ， 当 a = b 时 是 聚 光 情 形 。 下 图 是 望远镜 
系统 telescope 成 像 示 意图 。 


mo ee I oO — ~ 


Figure 4. telescope 成 像 示 意图 
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3.2 ”周期 性 磁 聚 焦 结构 (lattice) 
对 于 周期 为 L 的 周期 性 磁 聚 焦 结构 ， 假 定 Hill 方 程 为 
y" +RK(s)y=0,K(s+L)= K(s), (81) 


如 果 采 用 硬 边界 磁铁 条 件 ，K(s) 可 看 作 一 系列 阶梯 函数 的 组 合 。 在 储存 环 两 个 点 s, so 之 
间 ，Hil 方 程 的 解 y,y 可 写 为 


T ) = M(s/so) (4 i, | (82) 


其 中 
Mois) = (Gey sty) 6») 
以 及 
他 6 
连续 点 之 间 的 传输 矩阵 满足 级 联 方程 
M (82/89) = M (52/51) - M (51/50). (85) 
Xt JAPA LA ERR, HAEEREN 
M(s) = M(s + L/s), (86) 
个 周期 的 磁 聚 焦 结 构 ， 其 传输 矩阵 为 
M(s) = M(s+nL/s) = M(s)", (87) 


Min 一 00, WAS A Bi #4) 1 Fa KE BA REAR, ABA WR Re Fi EK). A A 
定性 条 件 ， 求 解 如 下 本 征 值 方程 


oy ae me 17 
mM (3 ) ea) (88) 
以 上 本 征 值 方程 有 非 零 解 的 条 件 是 


|M -A| =0. (89) 


DARRERE, HMMS = (7 h), ATTE AREE 
à? — (a + d)À + (ad — bc) = 0. (90) 
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上 式 为 一 元 二 次 方程 ， 其 中 ad — bc = |M| = 1， 它 的 两 个 解 入 ,和 2 满足 和 1 十 和 2 = a 十 
d, dy A2 = 1, 入 1， 入 2 可 写 为 

N12 = exp(tip) = cosu + isinp, (91) 
所 以 i i 

cosu = 5 (4 +d)= girace(M), (92) 


上 上 面 我 们 引入 了 一 个 新 的 实 变量 j, 称 之 为 相 移 变 量 ， 当 itrace(M) < 1 时 /为 实 变量 ， 
ts |trace(M)| < 1 为 稳定 性 条 件 。 在 /的 基础 上 ， 定 义 twiss 参 量 a, 8,7 如 下 


M = Icosp + Jsiny, (93) 
式 中 7 一 ( SF) r= ( 1) 为 单位 给 阵 。twiss 参 量 w, 9.1 的 表达 式 如 下 
a-d 
aa 2sinu’ (94) 
b = b/sinp, (95) 
y = —c/sinp. (96) 
由 |M| = 1 可 得 
By— oe =1. (97) 
和 矩阵 J 有 以 下 特性 
|J| = 0, trace(J) = 0, J? = —I. (98) 


对 于 任意 相 移 ji, wa 有 以 下 等 式 
(Icosuı + Jsinpy) - (Icos + Jsinplg) = Icos(py + u2) + Jsin(uı + u2) (99) 
那么 M 的 n 次 方 为 


M” = Icos(np) + Jsin(np). (100) 
总 之 ， 传 输 和 矩阵 M(s 可 由 四 个 参量 a(s), B(s), ys), 1 来 表示 
MO a (101) 
a(s), 6(s),Yy(5) 为 周期 函数 
a(s) = a(s + L), B(s) = B(s + L), y(s) = 7(s + L). (102) 


WISER. 
接 下 来 ， 讨 论 磁 聚 焦 结 构 a(s), 8(s),y(s) 的 传输 矩阵 。 可 以 证 明 yy? + 2ayy + By? = 
constant 为 不 变量 。 在 储存 环 上 两 个 点 s, so 处 有 等 式 


yy? + 2ayy’ + By” = your + 2aoyoyo + Boyd’, (103) 
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式 中 变量 y = y(s),y = y'(s), yo = y(S0), yb = y'(S0), 以 及 


Cr )= Ce sm) 
y Cc Ss yo J’ 
传输 矩阵 的 逆 和 矩阵 为 
i ( ‘o a 
所 以 
WASH, 


y0(S'y—Sy')? +2a0(S’y— Sy’) (—C'y+Cy’) + Bo(—C'y+ Cy’)? = yy? +2ayy' + By”, 


比较 上 式 两 边 的 系数 可 得 


B(s) = C? Bo = 25 Cao STe S?o, 
a(s) = —C C" Bo + (SC + SO — SS", 
(s) = C”? Bo — 29C" ao + S?o, 


p C? 25C S? Bo 
a |=| -CC SC+SC -SS ao |. 
y Cc”? —25'C' S52 Yo 


3.3 FODO 实 例 


用 和 矩阵 表达 为 


telescope 望 远 镜 系 统 是 没有 二 极 铁 的 FODO 磁 聚焦 结构 。 下 面 讨论 一 个 有 二 


(108) 


(112) 


极 铁 


的 FODO cell 磁 聚 焦 结 构 ， 忽 略 二 极 铁 的 弱 聚 焦 效 应 (二 极 铁 的 弯 转 半径 m])， 对 四 极 铁 


采用 莱 透镜 近似 ， 如 下 图 
F D 
| dipole i dipole 


e— 


L 
图 5. FODO 单元 
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I 
po Nl) 


四 极 磁 铁 的 焦距 为 上 FEFODO 单元 磁 聚 焦 结 构 的 周期 长 度 为 L。 由 和 矩阵 分 析 方 法 可 
得 FODO 单 元 相 移 /为 ， 
1L 


cosu = 1 — ga (113) 
O 在 聚焦 四 极 铁 中 间 最 大 ， 记 为 paz, 在 散 焦 四 极 铁 中 间 最 小 ， 记 为 Bwin, 其 表达 式 为 
nt 
Bmar — (114) 
sinu 
— sin” 
Bmin = Pm ma (115) 
sinu 


同样 ， 色 散 函 数 DD 在 聚焦 四 极 铁 中 间 最 大 ， 记 为 Daz, 在 散 焦 四 极 铁 中 间 最 小 ， 记 


iP ie sine 


7 三 SS x ) 116 
4p sin? (119) 
L? 1-— isin 


在 设计 中 ， 常 用 的 lattice 匹 配 程序 有 MAD、Transport 等 ， 以 下 是 一 个 MADX 匹 配 
的 FODO 单 元 。 


二 一 一 [ l [一 
21) FODO CELL EXAMPLE, 2022 MAD-X 5.08.01 22/11/22 08.30.46 
18. 


16. 


B (m). B (m). D (m) 


图 6. FODO lattice twiss 参数 图 
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4 AME 


带电 粒子 的 纵向 运动 讨论 粒子 能 量 加 速 的 问题 。 运 用 电场 加 速 带 电 粒 子 是 最 主要 的 
加 速 方法 ， 而 在 感应 加 速 器 、Betatron 加 速 器 中 运用 变化 的 磁场 产生 的 电场 来 加 速 用 的 
oe a 
ZARE, Te URE FT Pe AY ve 2] GR HS TD Bo BAT FET VO BAS TOPS 0 
Risk» CEL IY pet UU TY GR DR BER FT DN a VE Aa Ve FS FS A TT 1, 


图 7. BRAY AS ree ST] BR 


如 上 图 ， 假 定 高 频 间 际 宽度 g, 高 频 电 压 Vcoswt, 那么 高 频 电场 为 


V 
E, = —coswt. (118) 
g 


如 果 带 电 粒 子 的 速度 为 w 那么 粒子 通过 高 频 间 隙 时 获得 的 能 量 增益 为 


9/2 
Ap= f 1 (119) 
~g/2 g U 
sind /2 
eV JE eVT, (120) 


式 中 6 = w2 为 渡 越 角 ， T = PAURE, à 渡 越 因子 反应 了 加 速效 率 ， 一 般 为 了 
提高 加 速效 率 ， 可 降低 间 隐 长 度 9, 但 间隙 过 短 会 引起 打 火 ， 通 常 渡 越 因 子 值 约 在 0.8 附 
近 。 

加 速 间隙 是 加 速 谐振 腔 的 简化 ， 下 图 为 钟 形 单 cell、 双 cel 高 频 腔 的 示意 图 。 我 们 引 
入 分 路 阻抗 及 的 概念 


R,=— (121) 


Pare WA FEE tet SOTA] BR BES PS RIK RE PAA = R/L, LII BR w 
度 。 
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VS wien 
图 8. 典型 的 高 频 腔 示 意图 


4.1 ” 纵 问 运动 变量 的 选取 


粒子 纵向 运动 变量 指 单 粒子 动力 学 的 运动 变量 ， 纵 向 运动 本 质 上 当然 是 牛顿 运动 
方程 ， 它 是 二 阶 微分 方程 ， 我 们 选取 的 变量 适合 0 
把 纵向 运动 的 牛顿 方程 化 为 两 个 一 阶 微分 方程 ， 然 后 运用 哈密 顿 力学 方法 处 理 纵向 运 
动 。 


影响 环形 加 速 器 粒子 纵向 运动 的 因素 主要 有 高 频 腔 的 加 速 电 压 、 磁 铁 的 色散 效应 。 
我 们 主要 关心 的 是 粒子 一 圈 (one turn) 运 动 的 纵向 动力 学 ,也 就 是 研究 one turnia zA a] 
运动 变量 的 变化 ， 但 是 把 纵向 运动 变量 连续 化 、 平 滑 化 。 其 实 最 适合 研究 环形 加 速 器 
纵向 运动 的 方法 是 离散 化 的 一 圈 映 射 (one turn map)， 也 是 一 个 矩阵 化 的 分 析 方 法 。 对 
于 纵向 运动 ， 选 取 连 续 型 运动 变量 为 


AE 
fat (122) 


式 中 wo 是 纵向 运动 角 频 率 ，j 是 高 频 腔 射频 频率 的 谐 波 数 ，A 厂 为 经 过 高 频 腔 后 粒子 能 
量 的 增益 ，9, 为 同步 相位 ，9 一 内 为 每 次 粒子 经 过 高 频 腔 与 同步 相位 的 相位 差 。 


4.2 ”粒子 能 量 的 时 间 变 化 率 


假定 高 频 腔 的 射频 频率 了 = Aso, 高 频 电 压 幅 度 为 V, 粒子 的 电量 为 q, 那么 每 一 圈 粒 子 
经 过 高 频 腔 的 能 量变 化 为 


AE = qV sind, (123) 


从 而 粒子 能 量变 化 率 为 


~ AE- f = qVsind- f, (124) 


一 步 参考 同步 粒子 ， 可 得 以 下 一 阶 微分 方程 
d Eat oV 
dt^ hwo 2r 


在 上 式 中 ， 不 能 简单 将 wo 视 为 常数 ， 虽 然 我 们 将 高 频频 率 视 为 常数 ，wo 要 随时 间 变 
化 。 


(sing 一 sing,). (125) 
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4.3 ”粒子 运动 相位 的 时 间 变 化 率 


以 下 考虑 粒子 运动 相位 变化 ， 它 是 粒子 


为 7, 其 变化 率 为 


ab =. 


Be 里 . 


变化 的 函数 。 假 定 粒子 一 圈 运 动 时 间 


dr dC dv 
T OC v’ (126) 
式 中 ，C 为 粒子 运动 的 轨道 周 长 ，" 为 粒子 速度 。 轨 道 周 长 的 变化 是 由 于 磁铁 的 色散 效 
应 
dC dp 
C 中 a (127) 
ao 是 储存 环 的 动量 压缩 因子 ，z 为 粒子 动量 。 由 相对 论 运 动 学 有 以 下 等 式 
dp 1 dy 
pt ay eect 128 
p Paq’ oe 
dB 1 d 
式 中 6 为 归 一 化 速度 ，> 相 对 论 因 子 。 从 而 有 
dr dp db i d 
T mr = B = gz (Op z2) (130) 
dr ndy 
m (131) 
式 中 = oa —a, 为 滑 相 因 子 。 粒 子 相 位 变化 与 一 圈 运 动 时 间 的 关系 为 
dp dr 
z ee (132) 
do = ann, (133) 
从 而 
dọ -omn (134) 
ae dé AE 
= = db: fo = -hwo z huns p (135) 
最 终 我 们 得 到 粒子 运动 相位 的 时 间 变 化 率 为 
d = n 
PAG — bs) = — hwo Fag E 一 卫 。) (136) 
化 为 归 一 化 变量 可 得 最 终 的 方程 
d = 23 N) ,AE 
au? — bs) =—h “0 7B haa” (137) 


4.4 ”粒子 运动 的 哈密 顿 方程 
由 粒子 能 量变 化 的 时 间 变 化 率 、 粒 子 运 动 相位 的 时 间 变 化 率 得 到 了 两 个 一 阶 微分 方 


程 


d AE, qv 


d _ 72,2 0 (AE 
ge os) = ub eG) (139) 


以 上 两 个 方程 当然 可 以 化 为 一 个 二 阶 微分 方程 ， 可 以 得 到 一 个 单 摆 运动 方程 。 不 过 ， 

我 们 可 以 把 它们 化 为 哈密 顿 函数 、 哈 密 顿 方程 ， 运 用 已 有 的 哈密 顿 力学 分 析 方 法 来 进 

行 分 析 ， 比 如 相 空 间 体积 守恒 的 Liouville 定 理 ， 运 动 不 变 量 〈 如 能 量 守恒 等 ) 。 
EAE, 6 — 办 的 一 阶 微分 方程 ， 可 以 定义 以 下 哈密 顿 函数 


laa n AE WV a E, ; 
H = z” wo BB Tox? och! cosh 一 sinp, + cosh, + d,singds), (140) 
那么 哈密 顿 正则 方程 为 
d ,AE OH 
Par zp’ (141) 
d OH 
gO =~ 9088 (142) 


运用 以 上 哈密 顿 正则 方程 可 以 研究 粒子 在 纵向 相 空 间 运动 ， 存 在 基于 稳定 点 的 纵向 捕 
获 区 。 可 以 仔细 研究 纵向 捕获 区 的 大 小 ， 纵 向 捕获 区 的 分 界线 。 也 存在 基于 鞍点 的 不 
稳定 区 ， 在 相 空 间 不 稳定 区 ， 粒 子 作 发 散 运动 等 。 


5 讨论 


本 文 主要 讨论 了 环形 加 速 器 中 基本 的 横向 运动 和 纵向 运动 ， 重 点 放 在 了 横向 运动 。 
环形 加 速 器 最 基本 的 内 容 应 该 是 周期 性 磁 聚 焦 结 构 ， 磁 聚焦 结构 的 种 类 很 多 。 对 撞 
机 储存 坏 一 般 采 用 FODO 结 构 ， 同 步 辐射 光源 储存 环 一 般 采 用 MBA 结 构 ， 阻 尼 环 采 
用 TME 结 构 ，MAD 程 序 可 以 匹配 磁 队 焦 结 构 。 而 RLA/ERL 返 回环 需要 匹配 到 三 阶 束 
流光 学 ，Transport 程 序 可 以 匹配 到 三 阶 束 流 光学 。 环 形 加 速 器 粒子 动力 学 中 还 存在 很 
多 的 共振 问题 ， 横 向 运动 的 两 个 方向 之 间 ， 以 及 横向 运动 和 纵向 运动 之 间 存 在 耦合 。 
另外 ， 环 形 加 速 器 的 动力 学 孔径 也 是 一 大 类 问题 。 


References 


[1] S.Y.Lee, Accelerator Physics. 


20 


19 


20 


M.Reiser, Theory and Design of charged particle beam. 
周 衍 柏 ， 理 论 力学 教程 ， 高 等 教育 出 版 社 . 
郭 硕 鸿 ， 电 动力 学 ， 高 等 教育 出 版 社 。 


H.Goldstein, Classical Mechanics, Addison Wesley Publishing Company, Inc. 


J.D. Jackson, Classical Electrodynamics, John Wiley & Sons, Inc. 
金玉 明 ， 电 子 储存 环 物理 ， 中 国 科学 技术 大 学 出 版 社 。 

刘 乃 录 ， 加 速 器 理论 ， 清 华 大 学 出 版 社 。 

赵 籍 九 ， 尹 兆 升 ， 粒 子 加 速 器 技术 ， 高 等 教育 出 版 社 。 


Fifth general accelerator physics course, CERN Accelerator School, 1994. 
Fifth advanced accelerator physics course, CERN Accelerator School, 1995. 


A.Chao, M.Tigner, Handbook of Accelerator Physics and Engineering, World Scien- 
tific. 


S.Peggs,T.Satogata, Introducton to Accelerator Dynamics, Cambridge University 
Press & Peking University Press. 


朱 雄 伟 等 ，BXERL-FEL 及 其 试验 装置 -X-ray 一 机 两 用 的 探索 ，talk at IHEP, 2008. 


ARE, Design studies on the FEL-ERL at IHEP, talk at SLAC, 2009. 


AEB“, Design studies on the FEL-ERL at IHEP, talk at CESR, 2009. 


朱 雄 伟 等 ，Study on the characteristics of linac based THz light source, Chinese 
Physics C,33,10, 922 (2009). 


REIRE, Design study on the merger for BKERL-FEL, Chinese Physics C, 33, 8, 
696 (2009). 


H.Grote, et al., MAD-Methodical Accelerator Design, 2015. 


K.L.Brown, et al., Transport: A Computer Program for Designing Charged Particle 
Beam ‘Transport systems, 1983. 


Beam Dynamics of Circular Accelerator 


Zhu Xiong-Wei 
Institue of High Energy Physics, Chinese Academy of Sciences 


Abstract 


21 


In this paper, we present a pedagogical description of the single particle dynamics 
for the circular accelerator. The beam dynamics of circular accelerator can be 
classified into the transverse motion and the longitudinal motion. The materials of 
this paper come from the practice of the circular accelerator’s theory and design. We 
try to combine the physical problem and the mathematical method effectively and 
wholly. On the one hand, The requirement from the physical problem accelerate the 
developement of mathematics; On the other hand, more advanced the mathematical 
method, more easy the physical problem. One of the advantages is that we illustrate 
the beam dynamics of circular accelerator by Hamitonian Mechanics. 

Keywords: Circular accelerator, transverse motion, longitudinal motion, Hami- 
tonian Mechanics 
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